Strain localization in a shear transformation zone model for amorphous solids 
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We model a sheared disordered solid using the theory of Shear Transformation Zones (STZs). In 
this mean-field continuum model the density of zones is governed by an effective temperature that 
approaches a steady state value as energy is dissipated. We compare the STZ model to simulations by 
Shi, et al.\Y. Shi et al. PRL 98, 185505 (2007)], finding that the model generates solutions that fit the 
data, exhibit strain localization, and capture important features of the localization process. We show 
that perturbations to the effective temperature grow due to an instability in the transient dynamics, 
but unstable systems do not always develop shear bands. Nonlinear energy dissipation processes 
interact with perturbation growth to determine whether a material exhibits strain localization. By 
estimating the effects of these interactions, we derive a criterion that determines which materials 
exhibit shear bands based on the initial conditions alone. We also show that the shear band width 
is not set by an inherent diffusion length scale but instead by a dynamical scale that depends on 
the imposed strain rate. 
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I. INTRODUCTION 

Amorphous materials often comprise or lubricate 
sheared material interfaces and require more complicated 
constitutive equations than simple fluids or crystalline 
solids. They flow like a fluid under large stresses, creep or 
remain stationary under smaller stresses, and have com- 
plex, history-dependent behavior. Bulk metallic glasses, 
granular materials, and bubble rafts are just some of the 
disordered materials that exhibit these features [1, H, 0| ■ 

In this paper, we focus on strain localization, the spon- 
taneous development of coexisting flowing and stationary 
regions in a sheared material. Strain localization has 
been identified and studied experimentally in granular 
materials d, @ , bubble rafts [E 0] , complex fluids [1, Q , 
and bulk metallic glasses [lEf. Shear banding may play 
an important role in the failure modes of structural ma- 
terials and earthquake faults. Localization is a precursor 
to fracture in bulk metallic glasses [H and has been cited 
as a mechanism for material weakening in granular fault 
gouge on faults [111 ]. 

We model sheared material interfaces using the theory 
of Shear Transformation Zones (STZs). Spaepen first 
postulated that plastic flow could be modeled by track- 
ing localized zones or regions of disorder ■ Argon and 
Bulatov developed mean-field equations for these local- 
ized regions fl3l fl4| . Building on this work, Falk and 
Langer [l5[ introduced a mean-field theory which postu- 
lates that the zones, which correspond to particle config- 
urations that are more susceptible to plastic rearrange- 
ments than the surrounding particles, could reversibly 
switch between two orientational states in response to 
stress. Additionally, these zones are created and annihi- 
lated as energy is dissipated in the system. This theory 
captures many of the features seen in experiments, such 
as work-hardening and yield stress [ll| [TJ [ll| . 
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There are several complementary approaches to un- 
derstanding disordered solids. On the smallest scales, 
molecular dynamics and quasi-static simulations gener- 
ate a wealth of information about particle interactions 
and emergent macroscopic behavior — including shear 
banding P, |2(| — but they are limited to smaller 
numbers of particles and time scales. On the largest 
scales, phenomenological models such as viscoelasticity 
and the Dieterich-Ruina friction law [2l|, [52] , which has 
been studied extensively in the context of rock mechan- 
ics, describe stress-strain step and frequency responses, 
but to date these laws have not been derived from micro- 
scopic dynamics. 

A third approach focuses on the dynamics of collec- 
tions of particles and their configurations. Examples of 
this ap pro ach include the Spot Model for granular ma- 
terials IH, Soft Glassy Rheology (SGR) [23], and the 
theory of Shear Transformation Zones. The Spot Model 
generalizes dislocation dynamics in crystalline solids by 
postulating that a spot, a region of extra free volume 
shared among a collection of particles, executes a random 
walk through the material. Both SGR and STZ theories 
postulate that there are meso-scale configurational soft 
spots which are more susceptible than the surrounding 
material to yield under shear stress. 

Recently, Langer postulated that an effective temper- 
ature should govern the density of shear transformation 
zones (25|. Mehta and Edwards [2(| were perhaps the 
first to point out that although thermal temperature does 
not determine statistical distributions for macroscopic 
particles such as powders, the statistical properties of 
these systems could still be characterized by a small num- 
ber of macroscopic state variables, such as free volume. 
Cugliandolo, Kurchan, and Peliti [27J show that heat flow 
in these slowly stirred systems is determined by an effec- 
tive temperature. Several studies have confirmed that 
effective temperature can be used to characterize slowly 
sheared granular packings and bubble rafts [H, [29| . 

For many systems thermal temperature is not suffi- 
cient to cause configurational rearrangements, but slow 
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shearing causes the particles to ergodically explore con- 
figuration space. Although we make no attempt to calcu- 
late the entropy from first principles, we postulate that 
because the material is being sheared slowly, configura- 
tional (as opposed to kinetic) degrees of freedom are the 
dominant contributions to the entropy. The statistical 
distribution of steady state configurations should maxi- 
mize this entropy and be described by an effective tem- 
perature. 

The effective temperature is an important macroscopic 
state variable in this STZ description. It governs the 
properties of statistical distributions for configurational 
degrees of freedom in the material, such as density fluc- 
tuations. Langer hypothesized that shear transformation 
zones are unlikely, high-energy density fluctuations, and 
therefore the number of such zones in a system should 
be proportional to a Boltzmann factor which here is a 
function of the effective temperature instead of the usual 
thermal temperature [25| . In this manner, the density of 
shear transformation zones is related to the effective tem- 
perature and therefore the plastic flow is coupled to the 
local disorder. This is a mechanism for strain localization 
because a region with higher disorder is more suscepti- 
ble to flow under stress and flowing regions become more 
disordered. 

The effective temperature mechanism for strain local- 
ization is different from other proposed mechanisms gov- 
erned by thermal temperature or free volume. Griggs 
and Baker [3(| proposed that the strong temperature 
dependence of viscosity near the glass transition leads 
to thermal softening and strain localization. In that 
formulation the heat generated by plastic work raises 
the thermal/kinetic temperature, instead of the con- 
figurational/effective temperature. Lewandowski and 
Greer [3l| have shown that in a bulk metallic glass the 
thermal temperature diffuses too quickly to control shear 
localization in an adiabatic description, though Braeck 
and Podladchikov [32| have shown that a non-adiabatic 
thermal theory can explain the data. 

Alternatively, several authors have suggested that the 
free volume governs strain localization [lO, [33| . In these 
descriptions, plastic deformation in the solid is accom- 
panied by a larger free volume, which in turn softens 
the material. However, the exact relationship between 
the free volume and plastic deformation is difficult to 
determine. In this paper we suggest an alternative de- 
scription: the heat generated by plastic deformation is 
dissipated in the configurational degrees of freedom, ac- 
cording to the first law of thermodynamics. The local free 
volume is related to the effective temperature (for exam- 
ple, in a Lennard- Jones glass the particles can become 
more disordered by moving away from their equilibrium 
separations), but the effective temperature is the relevant 
macroscopic state variable governing particle dynamics. 

The STZ formulation in this paper is athermal — the 
thermal temperature T does not cause any configura- 
tional rearrangements. This is an approximation which 
is clearly appropriate for bubble rafts and granular mate- 



rials and may be relevant for metallic glasses well below 
the glass transition temperature. The athermal approx- 
imation significantly simplifies the STZ equations and 
clarifies the nature of the instability that leads to local- 
ization, but the resulting description does not capture 
features such as thermally activated creep or relaxation. 
As we discuss in Section TlVDi some behaviors of real ma- 
terials and simulations are likely due to thermal activa- 
tion or relaxation and will not be accurately represented 
in this athermal STZ formulation. 

In this paper we show how athermal STZ theory with 
effective temperature generates solutions with strain lo- 
calization. This paper is divided into sections as fol- 
lows: Section [IT] develops the constitutive equations for 
athermal STZ theory with effective temperature in a sim- 
ple shear geometry. Section IIIII compares numerical re- 
sults from the continuum STZ theory to simulations of 
Lennard-Jones glasses by Shi, et al. We compare macro- 
scopic stress-strain curves for different initial conditions 
and compare the strain rate and potential energy as a 
function of position and time inside the sheared mate- 
rial. Section IIVI investigates the stability properties of 
the STZ equations, and presents a description of how 
shear bands form and remain intact for long times. We 
show that an instability in the transient dynamics ampli- 
fies perturbations to the effective temperature, and that 
this process interacts with energy dissipation processes 
to determine the long-time behavior of the strain rate 
field. Section [V] concludes the paper with a discussion of 
our results and future directions. 



II. REVIEW OF STZ EQUATIONS 

A mean field theory for shear transformation zones, 
first derived b y Argon 1 131 , has been developed in a series 
of papers dHHH, OS |25| |. For completeness, this section 
reviews the athermal STZ equations and emphasizes their 
similarity to other models for elastic and plastic flow. 
Section PlI Al introduces and motivates the STZ equations. 
Section ITTBl tailors the equations to describe a particular 
geometry - an infinite strip of material driven from the 
boundaries at a velocity Vq ■ 



A. General equations 

When describing systems that undergo plastic and 
elastic deformation, it is convenient to break up the total 
stress tensor into hydrostatic and deviatoric components, 
(Ty = — pSij + ~pSij , because hydrostatic stress generally 
does not cause yield. To simplify notation, the deviatoric 
stress has been nondimensionalized by an effective shear 
modulus JI that characterizes the stiffness of the STZs. 
In STZ theory, the stress scale JI is important because it 
characterizes the stress at which the material begins to 
yield. In the slowly sheared materials we are modeling, 
the speed of sound in the material is very fast compared 
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to the rate of plastic deformation. In this case the stress 
gradients equilibrate very quickly, and we take the zero 
density limit of the momentum conservation equations. 
This results in static elastic equations for the stress: 



dai. 



= 0. 



(1) 



We further assume that the rate of deformation tensor is 
the sum of elastic and plastic parts: 
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where V/Vt is the material or corotational derivative de- 
fined as: 
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+ Au~uJkj — LJikAkj, (3) 



and Wjj = \/2{dvi/dxj — dvj/dxi). We also use the no- 
tation Aij — dAij/dt. 

The plastic rate of deformation tensor can be written 
in terms of dynamical variables from STZ theory. We 
postulate that under shear stress, each STZ deforms to 
accommodate a certain amount of shear strain, and can- 
not deform further in the same direction. This is modeled 
by requiring that each STZ be in one of two states: ori- 
ented along the principle stress axis in the direction of 
applied shear, which we will denote "+" , or in the per- 
pendicular direction, "— " . 

Under applied strain, the STZ will flip in the direction 
of strain, from "— " to "+" . Under shear stress in the 
opposite direction, the STZs can revert to their original 
configurations, which corresponds to a flip from "+" to 
"— " . These rearrangements or flips occur at a rate which 
depends on the stress R(s) and a characteristic attempt 
frequency, To. Because each STZ can flip at most once 
in the direction of applied strain, STZs must be created 
and annihilated to sustain plastic flow. Based on these 
considerations, the number density of STZs in each di- 
rection, n± obeys the following differential equation 



r n± = R(±s)n T -R(Ts)n ± +T (^e" 1 



-i/x _ 



(4) 



where R(±s)/to is the rate of switching per STZ as a 
function of stress, T is the rate at which energy is dissi- 
pated per STZ, and n^e^ 1 ^ is the steady state density 
of STZs in equilibrium. 

The first two terms in Eq. ((4]) correspond to STZs 
switching from "+" to "— " states and vice-versa, while 
the last term shows that the STZs are created at a rate 
proportional to n^e" 1 ^ and annihilated at a rate pro- 
portional to their density. The plastic rate of deformation 
tensor is given by the rate at which STZs flip: 
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Pechenik [l8[ generalized Eqs. ([J} and ([5]) to the case 
where the principal axes of the STZ orientation tensor 
are not aligned with principal axes of the stress tensor S{j . 
These generalized equations can be written in terms of 
two new variables. The first is A = n to t/n 00l where n tot is 
the tensorial generalization of and corresponds 

to the total density of zones in a sample. The second is 
/noo, where is the tensorial generalization 
of (n+ — n_) and corresponds to the STZ orientational 
bias. 

The rate at which STZs are created/destroyed, de- 
noted by r, is an important, positive quantity. Falk 
and Langer [l5| first proposed that T should be pro- 
portional to the rate of plastic work done on the sys- 
tem. Pechenik 18] refined this idea, noting that the rate 
of plastic work could be negative since energy could be 
stored in the plastic degrees of freedom. He proposed an 
alternative definition for T: the rate at which energy is 
dissipated per STZ, which is always positive. We adopt 
this important assumption in our STZ model. 

In previous papers, [HI, [2f|, the functional form of T 
was necessarily complicated because thermally activated 
switching can release energy stored in the plastic degrees 
of freedom. As noted by Bouchbinder et. al [l6j], the 
functional form of T is considerably simplified in an ather- 
mal description because an STZ can never flip in a direc- 
tion opposite to the direction of applied stress. Therefore 
no energy can be stored in the plastic degrees of freedom. 
The rate at which energy is dissipated per STZ is the rate 
at which plastic work is done on the system, Q, divided 
by the volume (area in 2D) density of STZs, eoA, where 
£o = Artec. In terms of these new variables, the rate of 
deformation tensor can be written as: 
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and the energy dissipated per STZ is 
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where s 2 = 1/2 s lj s tj , and C(s) = 1/2 (R(s) + R(-s)). 
Based on physical considerations, C(s) is symmetric func- 
tion of the stress that approaches zero as the stress ap- 
proaches zero and approaches a line with a slope of unity 
as the stress becomes large. Bouchbinder et al. [HI de- 
scribe a one-parameter family of functions with the cor- 
rect properties, and we will use the simplest of these: 
C(s) = -2+|s|+exp(-|s|)(2+|s|). 

Together, Eqs. ©,([7]), and the tensorial generaliza- 
tion of Eq. (TJ| describe the evolution of the density of 
STZs and the orientational bias of STZs. The density A 
is driven towards a Boltzmann distribution as energy is 
dissipated in the system: 



A = r 



(8) 
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The orientational bias of the STZs m,j is incremented 
to accommodate plastic strain, and reset as energy is 
dissipated in the system: 



T^ m ij _ 2 f,(—\r s ii 
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To close the system of equations we require an equation 
of motion for the effective temperature, x- As mentioned 
in the introduction, we postulate that the effective tem- 
perature describes the energy in the disordered configu- 
rational degrees of freedom. Ono et al. [28| show that at 
low strain rates the system is driven towards a specific 
steady state of disorder that corresponds to a specific 
value for the steady state effective temperature, Xoo- 

It is assumed that the effective temperature is propor- 
tional to the potential energy per particle, with a con- 
stant specific heat cq. This implies that the steady state 
effective temperature is also proportional to the energy 
dissipated by the configurational degrees of freedom in 
steady state. Energy balance requires that the effective 
temperature change at a rate proportional to the rate 
at which energy is dissipated in the system. However, 
this change in configurational energy must be limited as 
the effective temperature approaches its steady state. As 
X approaches Xoo, the heat generated by plastic work is 
dissipated by other mechanisms, which we do not track. 
Additionally, we assume that the configurational energy 
can be transferred between neighboring particles, so that 
there is a flux in the effective temperature which is pro- 
portional to its gradient. The resulting equation for ef- 
fective temperature is: 



these boundaries. The speed at the boundaries is Vq and 
the width is L. Therefore the average strain rate inside 
the material is 7 = Vq/L, while the local strain rate is 
denoted 7 = dv x /dy = 2D X °^ 1 . Assuming that there are 
no pressure gradients, Eq. ([2]) requires that the pressure 
field remain constant in time. 

The stress s and STZ bias m are traceless tensors. Let 
the off-diagonal elements of the tensors be denoted by 
s = s xy = s yx and m = m xy — m yx and the diagonal 



elements by sq = s a 



-s yy and mo 



By symmetry, all fields are constant in x, and Eq. (fT]) 
shows the stress tensor is constant in the y direction as 
well. Substituting this into Eq. ^ we find that the plas- 
tic rate of deformation tensor then has two independent 
components: 
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and the energy dissipation per STZ, T, is given by: 
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Because Sjj is independent of position, the tensor equa- 
tion of motion for the stress (Eq. J2])) can be integrated 
in y over the interval [y = —L, y = +L]. This results in 
simplified equations for s and sq: 



d 2 x 

coX = Q (Xoo - X) + D X~Qy~2' 



(10) 



where D x is a diffusion constant with units of area per 
unit time. 

In the athermal STZ formulation there is no mecha- 
nism for relaxation of the effective temperature to the 
thermal temperature T, so that the effective tempera- 
ture everywhere tends towards Xoo- This is in contrast 
to a thermal formulation that permits the effective tem- 
perature to relax in regions undergoing minimal plastic 
deformation. In Lennard-Jones and bulk metallic glasses 
well below the thermal glass temperature we expect any 
relaxation to be very slow, so that the athermal descrip- 
tion is a good approximation. However, even very slow 
relaxation processes may have a significant impact on lo- 
calization and we re-evaluate the athermal assumption in 
Section El 



B. STZ equations in an infinite strip geometry 

We will now derive the STZ equations of motion in 
the geometry of a two-dimensional infinite strip driven 
by the boundaries in simple shear. Let x be the direction 
of the velocity at the boundaries, and y be normal to 



M fV 2e s -\ Vo , 

so = -2^(^C(s)( S 4- mo )A)+s^, (15) 
fi \ tq s J L 

where (•) denotes an average over a field in the y- 
direction. For example, A is: 



K =kj_ L K{y)dy - 



(16) 



and D x ° y tal = Vq/{2L). Eq. © can be simplified to equa- 
tions for m and mo as follows: 

g-i/x 2 s 

m + m 7 = —Tm — 1 C(s)(- — m); (17) 

A T s 

e -!/x 2 s 

rn -m7 = -Tm — 1 C(s)(— -m ). (18) 

A t s 

In slowly sheared experiments and simulations, the strain 
rate 7 is always much smaller than the inherent at- 
tempt frequency I/to- Therefore the complicated ro- 
tation terms, m 7 and m/y, are very small and can be 
neglected. We rewrite the equations of motion so that 
times are in units of the inverse average strain rate, 
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1/7 = L/Vo, and lengths are in units of L. Note that 
stresses are already in units of the yield stress ~p. 
The resulting system of equations is given by: 

s = /i* (l - ^C(s) (J - m ) l\ - so; (19) 



S = -pL 



!!£ C (s)(££- mo )A) +.s: (20) 



-Tm- 



e^/x 2 



mo 



+ _C(3)( z -m)); (21) 
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where 



X = — reoA(xoo - X) + D x-E~2> 
c <V 



T = — C(s)(s — soTOo — sm). 
<7o 



(23) 
(24) 

(25) 



As noted in [l|| the density of STZs, eoA, is necessar- 
ily small. Eqs. ([19]), ([20]) and ([21]) each contain this 
factor in their numerators and they equilibrate very 
slowly compared to to, too and A, which are governed 
by Eqs. f2"T|). ([22]) and ([23]), respectively. Therefore we 
replace A, to and Too by their steady state values: 



A = A ss (x)=e-Vx ; 

, v / s/s, s < 1, 



(26) 
(27) 



too = mo(s, So ) = | SQ/ - 2j - >L (28) 

Below the yield stress the solid deforms only elastically 
because all the existing STZs are already flipped in the 
direction of stress. Three simple partial differential equa- 
tions remain: 
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(29) 
(30) 
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This set of equations has several dimensionless param- 
eters. The volume (area in 2D) fraction of the strip cov- 
ered by STZs or the probability of finding an STZ at a 
given point in the strip is eoA. A is the probability that a 
particle participates in an STZ. Assuming that STZs have 
a characteristic size of a few particle diameters, eo — An M 
is the characteristic number of particles in an STZ. The 



parameter fi* = /i//I is the ratio of average material stiff- 
ness to STZ stiffness, and corresponds to the slope of the 
linear elastic portion of a stress-strain curve. The specific 
heat in units of the STZ formation energy, en, determines 
how fast the effective temperature changes compared to 
plastic rearrangements. The parameter go = tqVo/L is 
the ratio of the two physical timescales in the system, the 
inverse STZ attempt frequency and the average inverse 
strain rate. Finally, Xoo is the steady state effective tem- 
perature in units of the STZ formation energy divided by 
the Boltzmann constant. 

The effective temperature diffusivity is the remaining 
parameter and requires some discussion. The dimension- 
less parameter D* in Eq. ([ITJ]) can be written as follows: 



D* x = (D x /c )(L/V )(l/L% 



(32) 



where D x — a 2 /td is a diffusivity with units of length 
scale squared per unit time. In STZ theory there are 
two candidates for the attempt frequency 1/tj> the lo- 
cal strain rate 7; oc and the inherent material attempt 
frequency 1/ro- The first possibility implies that the lo- 
cal potential energy diffuses at a rate proportional to the 
rate of local particle rearrangements, while the second 
implies that the energy flux depends on an inherent vi- 
brational/frictional timescale. 

Naively, one might think that the inherent frequency 
1/to is too high and results in a diffusion constant 
much too large to be physically reasonable. However, 
if one assumes that configurational energy diffusion is 
an activated process with large energy barriers (td = 
To exp(A-E/fcT)), then the resulting diffusion time td 
could be large and yet remain independent of strain rate. 

We have examined solutions to the STZ equations with 
an effective temperature diffusivity D* (a) proportional 
to the local strain rate and (b) independent of strain rate. 
In case (a) multiple shear bands are more likely, the ef- 
fective temperature takes much longer to reach its peak 
value, and the dynamics are more complicated. In case 
(b) a single shear band is more likely and the effective 
temperature quickly reaches its plateau. While further 
study is required to determine which diffusion frequency 
best describes the data, an effective temperature diffu- 
sivity that is independent of strain rate results in simpler 
dynamic equations, and we will use this assumption for 
the remainder of this paper. 



III. NUMERICAL INTEGRATION OF THE STZ 
EQUATIONS 

In this section we show that numerical integration of 
the simple, macroscopic STZ equations (Eqs. (|29H3ip ) 
qualitatively and quantitatively match molecular dynam- 
ics simulations of a Lennard- Jones glass by Shi, Katz, Li 
and Falk [l|. The simulations were performed in a sim- 
ple shear geometry with periodic Lees-Edwards boundary 
conditions. The parameter range for the simulations was 
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chosen to ensure that the glass exhibits strain localiza- 
tion, and stress-strain curves were computed for a variety 
of quenches and strain rates. 

The investigators calculated the potential energy per 
atom and strain rate as a function of position for several 
different strains. They used the data to investigate the 
hypothesis (found in the STZ model and other models) 
that the local plastic strain rate is related to an effective 
temperature or free volume. In [l|, they show that a 
value for the steady state effective temperature, Xoo, can 
be rigorously extracted from a thorough analysis of the 
simulation data. In this paper we numerically integrate 
the STZ equations in the same infinite strip geometry, 
and compare the results to simulation data. Building on 
the work of Shi, et al., we choose parameters for the STZ 
theory to match the conditions found in the simulations. 

In experiments and simulations, boundary conditions 
on the particles at the top and bottom of the strip im- 
pact the dynamics and may influence strain localization. 
Shi, et al. chose periodic boundary conditions for their 
simulations, while other investigators [lj| have chosen 
rigid rough walls. While both sets of simulations exhibit 
strain localization, features of the localization process are 
different in each case and it is unclear whether the dif- 
ferences are due to the dissimilar boundary conditions or 
other features of the simulations. 

Similarly, the STZ equations require boundary condi- 
tions for the effective temperature field at the top and 
bottom of the strip. We have studied the STZ equa- 
tions with both periodic and no conduction (dx/dy = 0) 
boundary conditions and found qualitatively and quanti- 
tatively similar results. For the remainder of this paper 
we will describe solutions to the STZ equations with pe- 
riodic boundary conditions for comparison to the results 
of Shi, et al. 



A. Matching parameter values to simulations 

The first task is to choose values for the parameters 
in STZ theory consistent with those in the simulations 
by Shi, et al. The simulation data for the stress-strain 
curves at a particular strain rate can be fit by a range 
of values for \oc, <7o, and the specific heat cq. However, 
a much narrower range of parameters fit both the stress- 
strain data and the data for the strain-rate as a function 
of position. 

We estimate an order of magnitude for each param- 
eter motivated by physical considerations. Table Q] lists 
the specific parameter values used in the numerical in- 
tegration. The steady state effective temperature \oo 
determines the steady state density of STZs: Aoo = 
exp[— 1/xoo]- To ensure that STZs are rare, should 
be around 0.1. Shi, et al., extracted Xoo = 0.15 from 
the data using the hypothesis that the strain rate is pro- 
portinal to e _1 / x , and we will use this value here. The 
parameter eo is the characteristic number of particles in a 
shear transformation zone. Studies of non-affine particle 



Parameter 


Value 


Xoo 


0.15 


eo 


10 


CO 


1 


D* 


0.01 


M* 


70 


<7o 


1 x 10 -6 


Xo 


0.68,0.69,0.74 



TABLE I: List of parameter values used in the numerical 
integration of the STZ model, Eqs. (|33l34jl . 



motion in a 2D Lennard- Jones glass [15j and a 2D system 
of hard spheres [34[ suggest that an STZ has a radius of 
a few particle diameters. This implies that the number 
of particles in an STZ is of order 10 for 2D simulations. 
Because cq is a dimensionless constant and does not de- 
pend on other scales in the problem, Co should be of order 
unity. The parameter /x* is the slope of the linear elastic 
part of the stress-strain curve plotted in units of the yield 
stress. It corresponds to the ratio of the elastic material 
stiffness to the STZ stiffness, and should be much greater 
than one. In the MD simulations, /i* is about 70. 

Eqs. (|29H31|) can be further simplified because (i* is 
large. For s < 1, the rate of plastic deformation is zero 
and and Eq. (|29| indicates that the stress increases pro- 
portionally to fi* . During this same time, Eq. (|30[) re- 
quires that So increases proportionally to s, which is less 
than unity. During the linear elastic response, the off- 
diagonal component of the stress is order fi* larger than 
the on-diagonal stress, which matches our physical intu- 
ition for simple shear. Numerical integration of Eqs. ((291 
13 1 1) confirms that sq remains about two orders of magni- 
tude smaller than s. We therefore use the approximation 
s = s and so = 0. This results in two simple equations 
for the STZ dynamics: 



( 1 - — C(s) (1 - m(s)) A 
(1-ms e /x Xoo-X 

co go 

d 2 x 



dy 



2 ' 



(33) 



(34) 



The remaining parameters are .D* and qo. As men- 
tioned in the previous section, we postulate that the ef- 
fective temperature diffusion D* is a constant indepen- 
dent of strain rate which we determine from the long-time 
diffusion of the shear band in the simulations. We find 
that D* ~ 0.01 matches the simulation data. 

The Lennard- Jones glass has a natural time scale 
to = &sl \/M/tsL, where osl is the equilibrium dis- 
tance between small and large particles, csl is the depth 
of the potential energy well between the two species and 
M is the mass. The average strain rate in the simulations 
varies from 2 x 10~ 5 in 1 to 5 x lO^ 4 ^ 1 . 



7 



If we assume that the STZ attempt frequency tq is 
approximately t , then q ~ 2 x 10~ 5 and the resulting 
numerical solutions to the STZ equations never exhibit 
localization. This is inconsistent with simulation data. In 
order to match the localization seen in the simulations, 
we are required to choose a value for tq that is an order 
of magnitude smaller than to- We find that qo ~ 1 x 
10~ 6 results in STZ solutions where strain rate matches 
results from the MD simulations. We will discuss the 
implications of this in Section IIV Dl 

Initial conditions for s and \ must also be determined. 
Because the stress is constant as a function of posi- 
tion and the simulations begin in an unstressed state, 
s(t = 0) = 0. The initial effective temperature x(Ui t = 0) 
is a function of position with many more degrees of free- 
dom. In the MD simulations, the initial potential energy 
is nearly constant as a function of position with fluc- 
tuations about the mean. Therefore, the initial effective 
temperature is also nearly constant with small amplitude 
perturbations about its mean. We expect fluctuations in 
the effective temperature to occur on the scale of a few 
particle radii in dimcnsionlcss units. The simulation box 
of Shi, et al. has a width of about 300 particle radii, so 
perturbations which span five particle radii have a nondi- 
mensionalized wavelength w = 1/60. 

The amplitude of the initial perturbations, 5xo, can 
be approximated empirically from the initial potential 
energy per atom as a function of position, shown in Fig- 
ured]^) . The potential energy at 0% strain exhibits small 
perturbations with standard deviation ~ 0.02 esL about 
a mean value of —2.51 esL, while the system reaches 
a maximum of —2.42 €sl in the shear band at larger 
strains. Assuming the maximum potential energy in the 
band corresponds to Xoo and the initial mean potential 
energy per atom corresponds to \o, the amplitude of ini- 
tial perturbations is about 8\ — 0.02. This is consis- 
tent with a thermodynamic calculation [3|| for the mag- 
nitude of fluctuations about the effective temperature, 

< (Sx) 2 >= rVM * o.oi. 

We use two different methods for generating an initial 
effective temperature distribution for the STZ equations. 
The first is to use a deterministic function with a single 
peak that serves as a nucleation point for shear bands: 
x{y,t = 0) = xo + Sxo sech(y/w). The second method 
generates a random number from a uniform distribution, 
and smooths those values using a simple moving average 
of width w. The resulting function is normalized so that 
its standard deviation from the mean value is <5xoi an d 
the mean is set at xo- 

The first type of initial condition is less physically re- 
alistic but more tractable because it generates at most 
a single shear band. The second type of initial condi- 
tion can generate solutions with varying numbers of shear 
bands, depending on the system parameters. We find 
that for the parameter range which best fits the simula- 
tion data, the STZ equations with random initial condi- 
tions generate a solution with a single shear band, which 
agrees with observations from the simulations and further 



validates our choice of STZ parameters. 



B. Comparison of macroscopic stress-strain 
behavior 

We numerically integrate Eqs. (|3"31 IM| with parameters 
values discussed above, for many different values of xo- 
FigureQJa) is a plot of stress vs. strain for three different 
initial preparations of a material starting from rest and 
driven in simple shear. Dashed lines correspond to MD 
simulation data from samples permitted to relax for three 
different amounts of time before being sheared. Solid 
lines are numerical solutions to the simplified STZ model 
(Eqs. (|3"3"|) and ([51)1 ) for three different average initial 
effective temperatures. 

All of the samples (for both the MD simulations and 
the STZ solutions) exhibit strain localization, which 
influences the stress-strain curves shown in Figure [T] 
We will discuss how localization affects the macroscopic 
stress-strain curves later in this section, and analyze the 
localization process in Section [TV] 

The STZ solutions shown differ only in the initial mean 
value for the effective temperature, xo, and are the best 
least-squares fit as a function of xo for each quench. Note 
that general features of the stress-strain curves match: 
there is a linear elastic segment, followed by a decrease 
in the stress slope as the material begins to deform plasti- 
cally, a peak at about 2 % strain, and stress relaxation as 
the material softens and stored elastic energy is released. 

The glass in the MD simulations behaves differently 
depending on how long the system was quenched before 
shearing, and the STZ model captures this behavior. A 
longer quench results in a more ordered solid and corre- 
sponds to a lower initial effective temperature. Numeri- 
cal integration of the STZ equations indicate that lower 
initial effective temperatures gcncrically result in higher 
peak stresses. This matches both the MD simulations 
and physical intuition: a more ordered solid takes more 
time and stress to plastically deform because more STZs 
must be created to permit the deformation. 

In Figure QJb) the slope (/x* ) of the linear elastic seg- 
ment of the most disordered glass (Quench III) is smaller 
than in the samples that were quenched for longer times 
before shearing. A linear fit to these data shows that 
/i* ~ 70 for Quenches I and II, while /i* ~ 60 for Quench 
III. 

We model the material surrounding the STZs as an 
elastic medium, and the parameter fi* is the ratio of the 
elastic material stiffness to the STZ stiffness. We assume 
that the STZ stiffness is constant for different sample 
preparations. Therefore a variation in /i* between sam- 
ples indicates that the surrounding clastic medium is less 
stiff in the more disordered material. This is consistent 
with the work of Maloney and Lemaitre, who have shown 
that the elastic shear modulus is smaller in more disor- 
dered materials due to non-affine particle motion (36j. 

Although we could better fit the MD simulation data 
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FIG. 1: (color online) (a)Simulation data and STZ theory 
for stress vs. strain from Ref. [J. Dashed lines are simu- 
lation data for three different initial preparations. The red 
curve corresponds to a sample quenched for the longest time 
(100,000 to)before the sample was sheared and has an ini- 
tial average potential energy per atom of -2.507. The green 
and purple curves were quenched for 50,000 to and 10,000 to 
respectively, and have initial PE/atom of -2.497 and -2.477. 
The solid lines correspond to STZ solutions with values for 
the initial effective temperature (xo) that best fit the simu- 
lation data: 0.062, 0.063 and 0.067, respectively. (b)Linear 
elastic regime of the simulation data for the three different 
quenches. 



by allowing fi* to vary between samples, in order to limit 
the adjustable parameters in the theory we fix fi* = 70 
for all samples. As expected, Figure QJa) shows that the 
numerical STZ results closely match the simulation data 
for Quenches I and II, while the best STZ fit systemati- 
cally deviates from the simulation data for Quench III. 



Comparison of strain localization inside the 
material 



over 800 % strain. This is in contrast to the much faster 
stress dynamics that attain steady state in less than 10% 
strain, as seen in Fig.[TJ In the numerical STZ solutions, 
the effective temperature attains a maximum in the same 
physical location within the strip as the strain rate. This 
is remarkably similar to the dynamics of the potential 
energy per atom in the simulations, and completely con- 
sistent with the assumption that \ is proportional to the 
potential energy. Shi, et al. were the first to systemat- 
ically check this hypothesis and they used it to extract 
various STZ parameters, such as Xoo [J- 




-1 -0.5 0.5 +1 

Nondimensionalized Y-Position (y/L) 



One feature of the infinite strip geometry is that it 
permits the system to achieve very large strains without 
fracturing, so that both theory and simulation can track 
system evolution over very large strains. The STZ solu- 
tions not only match the short-time macroscopic stress- 
strain behavior, but also match the long-time dynamics 
of strain localization within the strip. 

Figure [2] shows the strain rate as a function of position 
for various values of the strain (corresponding to different 
times) for (a) the MD simulations and (b) STZ theory. 
This figure also shows (c) the simulation potential en- 
ergy per atom and (d) theoretical effective temperature 
as functions of position. The simulation data is averaged 
over increments of 100 % strain, while the STZ model is 
evaluated for specific values of the strain corresponding 
to the midpoint of each binning range. 

Localization is evident in both the simulation data and 
numerical STZ solutions. The strain rate as a function of 
position inside the material exhibits very slow relaxation 



FIG. 2: (color online) (a) Simulation data for strain rate vs. 
position at various strains from Ref. [l[ . Strain rates are aver- 
aged over bins of 100 % strain. The dashed line corresponds 
to the imposed average strain rate, (b) Theoretical STZ solu- 
tions for strain rate vs. position at various strains. Strain rate 
is evaluated for specific values of the strain corresponding to 
the midpoint of each binning range for simulation data. The 
dashed line is the average strain rate, (c) Simulation data for 
potential energy as a function of position [lj] and (d) Theo- 
retical solution for the effective temperature as a function of 
position. The dashed lines correspond to the initial values of 
the potential energy per atom and effective temperature. 



D. Implications for constitutive laws 

Strain localization has a large effect on the short-time 
macroscopic stress-strain behavior of the system. For ex- 
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ample, Fig. [3] (a) shows two possible initial conditions for 
the initial effective temperature field. The average initial 
effective temperature is the same in both cases. One ini- 
tial condition for the effective temperature varies with 
position while the other is constant as a function of posi- 
tion. Figure[3](b) shows the resulting stress-strain curves 
for each case. When \o varies as a function of position, 
the system localizes and the resulting stress-strain curve 
does not reach as high a peak value, because the mate- 
rial releases elastic energy more quickly. The steady state 
stress is the same in each case. This is just one example 
of the general sensitivity of macroscopic state variables 
to microscopic details. 
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FIG. 3: (color online) Macroscopic stress-strain curves are 
sensitive to localization, (a) The initial effective tempera- 
ture distribution is either homogeneous (blue) or slightly per- 
turbed (red), though the average value xo is the same in both 
cases, (b) The stress vs. strain predictions are significantly 
different for these two different initial conditions. The per- 
turbed system localizes, as shown in Figure EJd), and this 
changes the resulting macroscopic behavior. 

The fact that localization influences the macroscopic 
stress-strain curves has important implications for con- 
stitutive laws. Constitutive laws provide a relation- 
ship between the strain rate, stress, and a set of state 
variables that characterize the internal structure of a 
material interface. They describe how features of mi- 
croscopic particle dynamics determine macroscopic fric- 
tional properties, and are used extensively in models of 
earthquakes [U HH, granular flows (37[, and machine 
control [38| . Several investigators have adapted STZ the- 
ory to describe the macroscopic behavior of lubricated 
surfaces 39|, dense granular flows [34[, and earthquake 
faults [40 1 , but all of these formulations assume that the 
STZ variables do not vary as a function of position across 
the sheared interface. Deriving new constitutive laws 



based on localized dynamics is beyond the scope of this 
paper. However, our analysis shows that localization is 
possible in STZ models and that the corresponding con- 
stitutive behavior is altered. Therefore it is important 
for localization to be included in modeling efforts, and 
this will be a direction of future research. 



IV. STABILITY ANALYSIS 

Continuum models for amorphous solids create a 
framework for understanding the instability that leads 
to localization in amorphous materials. In the previ- 
ous section we showed that STZ theory with effective 
temperature exhibits strain localization. In the following 
subsections we investigate the stability properties of the 
dynamical system given by Eqs. ([55)) and 



A. Steady state linear stability 

The first step in understanding how the simplified STZ 
model given by Eqs. ([3"3"[) and ([3"4"| responds to perturba- 
tions is to perform a linear stability analysis around the 
steady state solution. Setting s — in Eq. (|33[) . results 
in the following equation for the steady state flow stress 



1 



2eo 
9o 



C( S/ )(l-l/ S/ )exp[-l/ Xo 



(35) 



For the values of go and Xoo in Table [l] Sf = 1.005 - 
the yield stress is nearly equal to the steady state stress. 
The steady state solution for \ is spatially homogeneous: 
X(y) = Xoo- 

We now show that the perturbed system, s = Sf + 
s , x = Xoo + 5x(y), is linearly stable. To simplify nota- 
tion, we define the functions / and g: 



f(s, X ) - | 



9(s,x) = 



dx 
dt 



1 _ 2 fpC(£) (1 

qa 



to) 



2L 



dyj (36) 



^(l-m)e- 1/x ( X .- X ) 



go co 



V. 



dy 2 



(37) 



The perturbation to \ can be written as a sum of normal 
modes that satisfy the boundary conditions: 



$x{y) 



E 



Sxke 



riTT 

T' 



(38) 



The operator d 2 /dy 2 is diagonalized in the basis of nor- 
mal modes and therefore the dynamics of each normal 
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mode are decoupled. In the limit of infinitesimal pertur- 
bations, terms of order S\ 2 and s 2 can be neglected. This 
results in the following linear equations for each mode: 

df _ df 

s = -^{sf,Xoo)s + -^{s fl Xoo)Sxk, (39) 

SXk = -g^(sf,Xoo)s + -g^(s f ,Xoo)Sxk- (40) 

The second term in Eq. (|39p requires some additional 
explanation, because the operator df jds includes an in- 
tegral. The action of this operator on a normal mode 
$Xke iky is given by: 



df 
dx 



where 



= ^(s/,Xoo)^ J S Xk e m ^ L dy(41) 

(42) 



0, k ^0, 

F(sf, Xoo), k = 0, 



C{sj 



Sxk- (43) 



This analysis reveals a particularly interesting and im- 
portant feature of the STZ model dynamics which we will 
return to in the next section. There is a fundamental dif- 
ference between the dynamics of spatially homogeneous 
perturbations to x an d perturbations with zero mean. 
Because the stress is always spatially homogeneous, to 
linear order the stress dynamics depend only on the av- 
erage value of x and are completely unaffected by zero 
mean perturbations to x- 

In order to determine the linear stability of each mode, 
we calculate the two eigenvalues (as a function of k) for 
the steady state linear operator A ss (k), which is the two- 
by-two matrix given by: 



s 

&Xk 



— Ass 



A ss (k) 



s 

5Xk 



An(k) A 12 (k) 
A 2 i(k) A 22 (k) 



The two diagonal terms can be written: 
2 e - 1 /x=o £oAt * ( ( i 





s 




$Xk 



(44) 
(45) 



An 



1 C'( Sf ) 



C{s f ) 



A 22 = 



qa \\s f 
x 0( S/ - 1), (46) 
2eos/e~ 1//Xoc 
c go 

x (-C(s f )(l i)) G(s f 1) - D* x k\ (47) 

where is the unit step function. 

If either of the two eigenvalues of A ss has a positive 
real part for a particular value of k, then that mode grows 



exponentially and the system is unstable with respect to 
perturbations with that wavenumber k. 

For k ^£ 0, A\ 2 (k) is zero and the matrix is lower tri- 
angular, so the eigenvalues are simply An and A 22 . Be- 
cause C(s) is monotonically increasing and positive for 
s > 0, both eigenvalues are negative for all values of 
k ^ 0. The eigenvalues calculated for k = are negative 
also. Therefore, the STZ model is stable with respect to 
perturbations in steady state. An analysis by Foglia of 
the operator A 22 for an STZ model that included thermal 
effects resulted in the same conclusion 14111 . 



B. Time- varying stability analysis 

In simulations and numerical integration of the STZ 
equations, localization of strain first occurs before the 
stress reaches a steady state. Transient localization is 
also seen in numerical simulations of Spaepen's free vol- 
ume model by Huang et al. [12], and in the Johnson- 
Segalman model for complex fluids [43| . This motivates 
us to study the stability of transient STZ dynamics. 

The field x in the STZ model given by Eqs. (03 Dg]) can 
be rewritten as the sum of normal modes. As discussed 
in Section IIV Al for small perturbations with wavenum- 
ber k, the governing equations for the k = mode are 
fundamentally different from all other modes. This mo- 
tivates us to write the field x a t each point in time as the 
sum of a spatially homogeneous field x(t) and zero-mean 
perturbations: 



x(y,t)^x(t) + T, s ^( t y kv ^ 

k^O 



nn 

T (48) 



To analyze the transient stability of this system, we per- 
mit the k = mode, x(t), to be arbitrarily large but 
constrain the zero mean perturbations, Sxk, to be small. 
This is slightly different from the normal mode analysis 
in Section [TV A I where both k = and k ^ modes were 
small. Substituting Eq. (|4"5|) into Eq. (|36p and neglect- 
ing the second order terms in 5xk results in the following 
equation: 

k^O 

= /i* (l — 2—C(s)(l - m Q )e- 1 ^ 
V 1o 



x 1 + 



1 



1 



k^O 



6Xk{ty ky dy 



, e o, 



/U* ( l-2-C(s)(l-m )e" 1/x 
f(s,X)- 



(49) 



This indicates that to linear order, zero-mean perturba- 
tions to x do not affect the stress. Substituting Eq. 
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into Eq. (|3T[) results in the following linearized equation 
for each perturbation mode: 



dx 



{s(t),x(t),k) S Xk - 



(50) 



This is a linearized equation for the dynamics of small 
perturbations about a spatially homogeneous, time vary- 
ing solution. It is valid whenever the perturbations are 
small, even if the magnitude of x(t) is large. 

Physically, this corresponds to the following experi- 
ment: the system is started from a spatially homoge- 
neous initial condition and the system remains spatially 
homogeneous until a small, zero-mean perturbation in 
introduced at time r. The time-varying linear opera- 
tor §^(s(t),x(t)) describes the growth or decay of these 
small perturbations. The functions s(t), x(t) are the solu- 
tions to the STZ equations with spatially homogeneous 
initial conditions s(t = 0) = sq, x(y,t = 0) = xo- If 
the initial conditions are homogeneous, the STZ equa- 
tions are simply ordinary differential equations, which 
are much easier to solve than the inhomogeneous equa- 
tions. 

This description for the linearized start-up dynamics is 
relevant to many experiments and simulations. For ex- 
ample, the simulations by Shi, et al. begin from a state 
where the initial potential energy, which corresponds to 
the initial effective temperature, is roughly constant as 
a function of position. This corresponds to a spatially 
homogeneous initial condition for the effective tempera- 
ture, x{t = 0, y) = Xo = c. Additionally, the MD samples 
are started from an unstressed state, which corresponds 
to so = 0. The potential energy does vary slightly as 
a function of position, which corresponds to small, zero- 
mean perturbations to x introduced at time t = 0. 

We define the stability exponent uj c (k, t) in the follow- 
ing manner. Let s c {t),x c (t) be the unique spatially ho- 
mogeneous solution to the STZ equations starting from 
the initial condition sq = 0, xo ~ c. Then the exponent 
at time r is defined as: 

Qg 

u c (k,T) = ^(s c (r),x c (r),fc) 



2e- 1 /x« 



e s c 



Co 9o 



xl 



x6(s c -l) 



--1 C( Sc )(l--) 
D* x k\ 



(51) 



where the functions s c (t) and x c W are understood to be 
evaluated at t = r and Q(s) is the unit step function. 

This exponent describes the rate of growth or decay 
of a small perturbation with wavenumber k introduced 
at time r to the solution s c (t),x c (£)- If the real part of 
uj c (k,t) is greater than zero then the perturbations grow 
exponentially and the system is unstable, and otherwise 
the system is stable with respect to perturbations. 

Note that oj c (k,t) contains a single term that depends 
on the wavenumber k — this term is proportional to 



— k 2 and corresponds to the diffusion of effective tem- 
perature. Because diffusion can only act to stabilize per- 
turbations to x, the most dangerous mode corresponds 
to k — ±tt/L. 

A plot of uj c (k,t) for c = xo = 0.6 and k = tt/L is 
shown in Fig. 0] The system is marginally stable with re- 
spect to perturbations during the linear elastic response 
of the system, highly unstable at intermediate times, and 
becomes stable with respect to perturbations as x(t) ap- 
proaches Xoc- A lower bound for stable x is determined 
from Eq. (|5Tj) as follows. Because the diffusion term is al- 
ways negative, an upper bound on the stability exponent 
for all values of k is given by: 



LO c (k,t) < H(s)e~ 



-1/XcO) ( Xoo - Xc(t) 

xl(t)-i 



(52) 



where H(s) is zero for s < 1 and positive for s > 1. 
Therefore the stability exponent is less than zero when 
x{t) satisfies the following inequality: 



> 
which defines 



Xc 



x(t) 2 



- 1, 



(53) 
(54) 



X(t) > - (-1 + Vl + 4Xoo) =x crit .(55) 
This analysis indicates if the initial effective tempera- 
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FIG. 4: The stability exponent ui c (k,t) with c = 0.06 and 
k — tv/L as a function of total strain for spatial perturbations 
to the effective temperature in the limit k — > 0. The system 
is linearly unstable with respect to spatial perturbations to 
X when Lu c (k,t) is greater than zero, and stable otherwise. 
The figure contains two scales for the average strain in or- 
der to simultaneously illustrate both the short- and long-time 
behavior. 

ture is above a critical value, xo > Xcrit' t nen small 
amplitude perturbations are stable and the system will 
not localize. However, if the initial effective tempera- 
ture is less than this value (plus a diffusive term that de- 
pends on the wavenumber), there exists a period of time 
where small perturbations to the effective temperature 
will grow. A time-varying stability analysis of thermal 
STZ equations by Foglia shows that the effective tem- 
perature is transiently unstable to perturbations in that 
context as well [4l|. 
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This description of perturbations to a spatially homo- 
geneous, time-varying, "start-up" trajectory is remark- 
ably similar that given by Fielding and Olmsted [43| in a 
study of shear banding in the Johnson-Segalman model 
for complex fluids. The stability of the solution to per- 
turbations is simpler to analyze in the STZ formulation 
than in the complex fluid model because the linearized 
evolution of perturbations to x are decoupled from the 
stress, as shown in Eq. (|49|) . This suggests that the mech- 
anism for shear band nucleation in "start-up" flows may 
be similar in a variety of materials. 

C. Finite amplitude effects and a criterion for 
localization 

While linear stability analysis determines when pertur- 
bations grow, by itself it does not provide enough infor- 
mation to determine whether or not the integrated STZ 
equations exhibit strain localization. Many of the numer- 
ical STZ solutions for an unstable initial effective temper- 
ature do not localize if the initial perturbations are "too 
small" . For example, numerical integration of the STZ 
equations with an unstable initial value for the effective 
temperature xo = 0.09 < Xcrit — 0-13 and a small ini- 
tial perturbation Sxo = 0.001, result in a solution that 
is virtually indistinguishable from a homogeneous sys- 
tem. The purpose of this section is to derive a criterion 
based only on initial conditions (the mean value of the 
initial effective temperature xo and the amplitude of ini- 
tial perturbations Sxo) that determines which materials 
will exhibit long-time strain localization. 

At first glance, the fact that a finite perturbation is re- 
quired to generate localized solutions seems to contradict 
our assertion that the system is linearly unstable. How- 
ever, linear stability equations are only valid for short 
times when perturbations are small - they give no in- 
formation about the long-time behavior of the unstable 
perturbed/inhomogeneous states. Nonlinear system dy- 
namics can enhance or inhibit localization in the inho- 
mogeneous states. Upon further examination, we recog- 
nize that the "finite-amplitude effect" is due to nonlin- 
ear dynamics involving the interaction of two dynamical 
processes: perturbation growth and energy dissipation. 
Based on this understanding we derive a criterion for 
which initial conditions result in shear banding. 

We first show that rate of energy dissipation as a func- 
tion of position determines whether perturbed states re- 
main inhomogeneous or become uniform. Our guiding 
intuition is that perturbations that grow quickly permit 
most of the energy to be dissipated in the incipient shear 
band, which enhances localization. Otherwise the energy 
is dissipated throughout the material, which inhibits lo- 
calization. 

Eq. can be rewritten in terms of the rate of plastic 
work dissipated by the system, Q: 

a = 1 - -, (56) 



where Q is spatial average of the rate of plastic work 
Q(y). Q can be resolved into a product of two terms 
- the rate of plastic work per STZ, T(s), that depends 
only on the stress and is spatially invariant, and the den- 
sity of STZs, A(x), which depends only on the effective 
temperature and is a function of position: 

Q = —sC(s)(l - m)e- 1 ^ = r(i)AM. (57) 

In steady state the average Q = r(s)A is constrained to 
be a constant, but the value Q(y) — T(s) exp [— l/xiv)] 
is not similarly constrained. This admits the possibility 
that the plastic work dissipated is very small at some 
positions and large at others, and this must occur for the 
system to remain localized. 

To see this, first note that if the effective temperature 
is below Xoo , Eq. (fTOf requires that the effective temper- 
ature must always increase proportional to Q(xoc — x) 
or diffuse — there is no relaxation towards thermal tem- 
perature in this approximation. Therefore, a spatially 
heterogeneous effective temperature distribution will be 
nearly stationary only if Q is very close to zero whereever 
X 7^ Xoo- In other words, strain localization will only per- 
sist if the rate of plastic work dissipated outside the shear 
band is very small. 

The rate at which plastic work is dissipated is propor- 
tional to exp [— l/x(y)]- A rough assumption about the 
nonlinear, inhomogeneous effective temperature dynam- 
ics is that the largest initial perturbation of amplitude 
Sxo continues to grow at the rate we derived using linear 
stability analysis: u>(t)8x{t), while the effective temper- 
ature in regions outside the incipient shear band grows 
at the rate of a homogeneous system, x = / (x(*)j s (t))- 
This assumption is consistent with behavior seen in nu- 
merical solutions to the STZ equations. 

Figure [5] shows the behavior of the effective temper- 
ature as a function of time for (a) a system that local- 
izes and (b) a system that does not localize, illustrating 
the two different initial stages. In both plots the upper 
solid line is the maximum value of x(y) inside the shear 
band, the lower solid line is the minimum value of x(u) 
outside the shear band, and the dashed line is a homo- 
geneous solution with the same mean value, %. In plot 
(a), the perturbation grows quickly compared to the ho- 
mogeneous solution for x, while in (b) the perturbation 
grows slowly. 

We are interested in whether the system dynamics fo- 
cus energy dissipation inside the shear band or spread 
energy dissipation evenly throughout the material. Ide- 
ally, we would like to compare energy dissipation inside 
the band (Qi n ) to energy dissipation outside of the band 
(Qout) at various times t. However, we only have informa- 
tion about the initial conditions, xo an d <5xo- Therefore 
we compare the time derivatives of Q inside and outside 
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ty t max 



FIG. 5: (color online) Plot of effective temperature \ as a 
function of time for (a) an STZ solution for Sxo = 0.01 that 
localizes and (b) an STZ solution for Sxo = 0.001 that does 
not localize. The dashed blue line indicates the homogeneous 
solution for x as a function of time. The red(orange) solid 
lines corresponds to the function x(y) evaluated at a point 
located within (outside of) the shear band. The shaded region 
indicates the range of values spanned by x(y) a * every time, so 
that a large shaded region for long times indicates a long-lived 
localized state. Dashed vertical lines correspond to calculated 
values for the time the material first reaches the yield stress, 



the band: 



dQ i7 
dt 



dt 



T(s) 



e- 1/x >~ d X ir, 
(x») 2 dt 



r( s y 



= r( s )- 



;-l/X0 
,-1/Xo 



w X0 (» S x (t); 
dxout 



{Xout) 2 dt 

p-i/xo 

(xo) 2 



(58) 



(59) 



We approximate Xin and Xou t as the initial condition, X o- 
This approximation is valid only for short times when X 
has not changed appreciably from X o inside or outside 
the band. The ratio of the derivatives is given by: 



u X o(k,t)6 X (t) 



(60) 



Equation ((60|) indicates that for short times when X has 
not changed appreciably, the energy dissipation is deter- 
mined by comparing the growth rate of perturbations to 
the growth rate of the mean effective temperature. 
Equation (|60|) must be evaluated at a particular time 



and wavenumber. We evaluate u> Xo (k,t) for the most un- 
stable mode, which corresponds to k = ir/L. For nota- 
tional simplicity this value of k is assumed in the follow- 
ing derivation. In writing Eqs. (|58[) and (|59|) . we required 
that the time t is small enough that X remains close to 



Xo- Numerical solutions confirm that X changes very lit- 
tle as s increases from zero to its maximum shear stress. 
This stress, s m , can be approximated as the solution to 
Eq. (J37J) with s = and \ = Xo- 



= 1 - 2— C(a m )(l 
<?o 



1 



-l/xo 



(61) 



We must also approximate the numerator of Eq. [5D] 
at the time when the perturbations are growing most 
rapidly. Naively, one might expect this to be the maxi- 
mum of lo Xo (t) times the amplitude of the initial pertur- 
bation, 5 X o- This underestimates the numerator because 
the amplitude of the perturbation increases with time. 

A better approximation can be found by noting that 
in numerical solutions, u> Xo (t) rises almost linearly from 
zero at the time the material reaches the yield stress to 
a maximum at nearly the same time as the material at- 
tains its maximum stress, s m . This behavior is seen in 
Fig. [4j Equation (|51[) shows that lu Xo (t) depends on t 
only through s(i) and x(t), so we evaluate the stability 
exponent at (s m , Xo)- Therefore we can approximate the 
linearized equation for the perturbations, Eq. (|50[) . as: 



Sx(t) 



u Xo (s m ,xo) 



{t-ty)6 X {t), t m >t> ty, (62) 



where t y is the time when the material first reaches the 
yield stress and t m is the time at the maximum stress. 
The solution to this differential equation, evaluated at 



0X{t m ) = 5xo exp[ a —\. (63) 

Therefore a lower bound on the numerator is given by: 



max [u^ 



Xt)5x{t)\ > u(s m ,Xo)S X o exp 



w{s m ) At 



(64) 



where At = t m —t y is nearly constant in all the numerical 
STZ solutions and is about 0.03 in units of strain. 

While better than the first guess, this estimate is likely 
to be low because u Xo (t) is generally not sharply peaked 
about t m , so the maximum growth rate occurs at later 
times than t m . To account for this and to match the 
numerical solutions, we choose At = 0.10. 

Substituting these approximations into Eq. 1601 we de- 
fine localization ratio 1Z as follows: 



n$x,xo) 



1 



>e C(s m ) 



Vxo 



,(65) 



where w(s m ,xo) is defined by Eq. (|5T|) . 1Z is large if the 
system tends to focus energy dissipation inside the band 
and small otherwise. 

We now systematically compare the value of the local- 
ization ratio to the degree of localization in a numerical 
STZ solution, and find that the ratio is an excellent cri- 
terion for localization. 
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For each set of initial conditions, (xo, ^Xo)j we com- 
pute a numerical solution to the STZ equations and then 
calculate the Gini coefficient for the strain rate as a func- 
tion of position for each point in time. The Gini coeffi- 
cient, (j)(t), is a measure of localization and is defined as 
follows 144, 451: 



*w - ^EEi^')-^W) 



,2 A 1 



-i/xte,*) 



2n 2 A 



, (66) 



where the function x is evaluated at n points in space, 
{y n }- A discrete delta function has a Gini coefficient 
equal to unity, and a constant function has a Gini coeffi- 
cient equal to zero. We define the Numerical localization 
number $ as the maximum value of <p(t) over t. If $ is 
close to unity the strain rate is highly localized (as in 
Figure H{b)) an d when <I> is close to zero the strain rate 
remains homogeneous. 

Figure [B] shows that the numerical localization number 
$ (rectangles) and the theoretical ratio 1Z (contour plot) 
exhibit the same behavior as a function of xo and Sxo- 
This indicates that the theoretical 1Z is a good criterion 
for the numerically computed strain localization $. 

The shading of the rectangles corresponds to the value 
of <E> for a given set of initial conditions: red corresponds 
to highly localized strain rate distributions (<f> > 0.8), 
yellow to partially localized solutions (0.8 > $ > 0.3), 
and dark blue to homogeneous solutions ($ < 0.3). The 
contour plot corresponds to values of 7Z: light pink for 
1Z > 0.6, darker shades for 1Z < 0.6. The data in figure[6] 
show that to a good approximation, JZ > 0.6 corresponds 
to solutions with strain localization, and 1Z < 0.6 corre- 
sponds to homogeneous solutions. 

The localization ratio 1Z is a function only of the STZ 
parameters and the initial conditions for the function 
x(y)- This is potentially a very useful tool for incor- 
porating effects of localization into constitutive laws. It 
determines an important material property without re- 
quiring computation of the full inhomogeneous system 
dynamics. 



D. Long-time behavior 

The localization ratio compares the energy dissipated 
inside the shear band to the energy dissipated outside 
the shear band, and determines whether the nonlinear 
dynamics enhance or inhibit shear band formation. An 
interesting and more difficult question is what determines 
whether the system dynamics result in a single shear 
band or multiple shear bands. The MD simulations gen- 
erally have a single shear band except at the highest 
strain rates (see Figure 2 in [l[). In the parameter range 
used to match the MD simulations, the STZ solutions 
also exhibit a single shear band. However, outside this 
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FIG. 6: (color online) Comparison of numerical localization 
$ to theoretical 1Z as a function of xo and log 10 ((5%o)- Each 
rectangle in Figure [5] corresponds to a solution to the STZ 
equations for a particular choice of (xo, Sxo)- The shading of 
the rectangle corresponds to the value of $ for a numerical so- 
lution with the indicated initial conditions: red corresponds to 
highly localized strain rate distributions, 1 > $ > 0.8, yellow 
corresponds to partially localized solutions 0.8 > $ > 0.3 and 
dark blue corresponds to homogeneous solutions 0.3 > $ > 0. 
The background shading corresponds to a contour plot for 
1Z. Light pink corresponds to values of 1Z > 0.6, dark purple 
corresponds to values of TZ < 0.1, and contours are at 0.1 
intervals in between. The yellow rectangles consistently over- 
lay the yellow contour, showing that the theory and numerics 
agree on the localization transition. 



parameter range it is possible for the STZ equations to 
generate solutions with multiple shear bands. 

It appears that the STZ equations generate a single 
shear band when the system is in a parameter range 
where perturbations to x are highly unstable. In this 
case the largest amplitude perturbation grows rapidly, 
resulting in a large rate of energy dissipation at that po- 
sition. Consequently the energy dissipation rates outside 
the shear band are smaller, which inhibits the growth of 
smaller amplitude perturbations. A full analysis of these 
dynamics is beyond the scope of this paper. 

However, we use the fact that there is a single shear 
band to determine the width of that band. As discussed 
above, most of the strain must be accommodated in the 
shear band if the strain is to remain localized. Xoo sets 
the steady state of disorder in the system, and we assume 
that x attains this value in the shear band. We can 
estimate the width of the band as a function of time, 
WE{t), by postulating that all of the shear band strain is 
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accommodated plastically in the shear band: 

= i-^ L fiy pl {y)dv, (67) 

o = i- ^l C ( s (t))(i- _L) e -vx~ ;(68) 

2L g V SV-) J 

WE(t) go cxp[l/xoo] . . 

2L 2C(s(t))(l-^Y l ° yJ 

where s(t) is the numerical solution to the STZ equation 
for the stress as a function of time. 

We also evaluate the width of the shear band in the nu- 
merical STZ solution for the strain rate. The numerical 
shear band width (wN(t)) is computed as the width of the 
strain rate vs. position curve at a value corresponding to 
the imposed strain rate, 1 = Vq/L. Figure \7\ compares 
the numerical shear band width to the estimated width, 
as a function of time. The two are in agreement. The nu- 
merically computed width is systematically larger than 
the estimated width because we chose a relatively low 
cutoff, 1 = Vq/L, for the strain rate in the band. 

The estimated width is not a prediction for the width 
of the shear band because it depends on the numerical 
solution for the stress as a function of time. However, it 
does suggest that because the STZ equations contain no 
natural length scale for the width of a shear band, the 
system dynamically sets the width based on the imposed 
strain rate. 

The shear band has a well-defined effective tempera- 
ture (xoo) and therefore accommodates a fixed rate of 
plastic strain per unit width at a given stress. Together, 
the stress and the imposed strain set the shear band 
width: the shear band must grow to the width required 
(at a fixed stress) to accommodate all of the imposed 
strain. While we cannot compute the shear band width 
analytically, we have recast the problem in terms of a 
potentially simpler one to solve: determining the stress 
relaxation as a function of time. 

The value of the steady state stress has a large effect 
on strain localization. While the spatial dependence of 
Q(y) is through x(y), the magnitude of Q depends on 
the value of the stress through T(s). When the system 
is in the flowing regime, the nondimensionalized stress s 
is always greater than one. The energy dissipation per 
STZ, r(s) oc sC(s)(l — 1/s), approaches zero as the stress 
approaches one from above. Therefore, if the steady state 
stress is very close to unity, the effective temperature 
dynamics become very slow - they are "frozen in" by 
the stress dynamics. In this parameter range the rate of 
plastic work can become very small at some positions as 
the system approaches steady state and this is exactly 
what permits a localized strain state to be long-lived. In 
contrast, when the steady state stress is too large the 
magnitude of Q is also large (regardless of x(u)) an d the 
material will not remain localized. 

Numerical integration of the STZ equations confirms 
that localization occurs only when the steady state stress 




Strain 

FIG. 7: (color online) Data for shear band width as a function 
of strain. The red dashed line is the numerical width wn 
calculated from the STZ solution for the strain rate. The solid 
green line is the estimated shear band width we calculated 
from the STZ solution for the stress only. 



is close to unity. STZ solutions exhibit localization when 
go is 0(1O~ 6 ) and Xoo is 0(1/10), which corresponds to 
a steady state stress that is very close to the yield stress 
(sf = 1.005). However, they do not exhibit long-lived 
shear bands if the system is driven more quickly (go ~ 
0(1O -5 )), which corresponds to a steady state stress of 
S/ ~l.l. 

This helps us interpret MD simulation data at higher- 
strain rates. In addition to the data shown in Fig- 
ure [TJ Shi, et al. have computed stress-strain curves 
for two higher strain strain rates (five and 25 times 
greater, respectively, than the strain rate used to gen- 
erate Fig.rjj m. 

The lowest strain rate is denoted VI and the higher 
strain rates are denoted V2 and V3. In this paper we 
defined the timescale ratio go using VI. Therefore V2 
corresponds to 5go and V3 corresponds to 25go- 

We numerically integrate the STZ equations at 5go and 
25go and find that the resulting stress-strain curves do 
not match the simulation data. The theoretical STZ 
steady state flow stress is the same for all three strain 
rates, while in the MD simulations the flow stress for a 
system driven at V3 is significantly (~ 15%) higher. 

One way to eliminate this discrepancy is to assume 
that we have misidentified the original timescale ratio go. 
For a given value of go (and Xoo), we can calculate the 
flow stress s / using Eq. (|35[) ■ To match the simulation 
data, the flow stress that corresponds to 25go should be 
15% higher than the flow stress that corresponds to go- 
In Table |TT] we have calculated the flow stress ratios for 
several values of go- This table shows that go ~ 1 x 10~ D 
better fits the MD simulation data for different strain 
rates than g ~ 1 x 10~ 6 . The higher value for g also 
makes sense physically because it corresponds to an STZ 
timescale tq which is simply the natural timescale for the 
Lennard- Jones glass, to. 

Our athermal model requires a small value for go to 
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go Sf(5qo)/s f (go) S/(25qo)/s/(qo) 

1.6 x 10~ 6 1.005 1.028 

1 x 10~ 5 1.028 1.144 

1 x 10" 4 1.185 1.695 



TABLE II: Comparison of the STZ steady-state flow stress 
for different values of go. 

generate solutions with long-lived localized states, but a 
larger value of qo to be consistent with MD simulations 
data for different strain rates. This disparity is likely 
an indication that the athermal approximation does not 
adequately capture long-time behavior of Lennard-Jones 
glasses. 

In an athermal description, the effective temperature 
always tends towards Xoo because there is no thermal 
relaxation. Alternatively, if the effective temperature 
can relax towards a thermal bath temperature in regions 
where there is little plastic deformation, localized states 
can persist even when st is large. It seems likely that 
a thermal description would generate localized solutions 
with higher values for qo- Although the athermal STZ 
formulation reproduces many aspects of the MD simula- 
tions, a model with thermal relaxation may be required 
to capture more details of the MD simulations. This will 
be a direction of future research. 



V. CONCLUDING REMARKS 

We have investigated the athermal STZ model with 
effective temperature for amorphous materials in a sim- 
ple shear geometry. In contrast to other studies, the 
effective temperature field varies as a function of posi- 
tion and by numerically integrating the STZ equations 
we have shown that the resulting solutions can exhibit 
strain localization. The numerical STZ solutions match 
the stress-strain curves for MD simulations of Lennard- 
Jones glasses by Shi, et al. and exhibit strain rate and 
effective temperature fields that are consistent with those 
seen in simulations. 

The continuum STZ model provides a framework for 
understanding how shear bands nucleate and persist for 
long times. We have shown that the model is unstable 
with respect to small perturbations of the effective tem- 
perature during the transient dynamics, though it is sta- 
ble to perturbations in steady state. Interestingly, shear 
localization does not always occur when the system is 
transiently unstable because shear bands can only form if 
the perturbations grow quickly compared to the homoge- 
neous solution. This is a nonlinear effect which is driven 
by uneven energy dissipation in the system. Using rough 
approximations to these nonlinear dynamics we were able 
to derive a localization criterion that depends only on the 
initial conditions for the effective temperature. 

In the parameter range studied, STZ theory predicts 
that the sheared material attains a state with a single 



shear band. The effective temperature field reaches its 
steady state value Xoo in the shear band and remains 
close to its initial value outside the band. To a good 
approximation, all of the strain imposed at the bound- 
aries is accommodated in the band, and we determine the 
width of the shear band using the STZ expression for the 
plastic strain rate when x = Xoo- This analysis implies 
the shear band width is not determined by an inherent 
length scale in the material, but instead by a dynamical 
scale set by the imposed strain rate. 

Interesting questions remain concerning the interaction 
between shear bands and material boundaries. We eval- 
uate the STZ equations with simple, periodic boundary 
conditions on the effective temperature field. In our nu- 
merical solutions, the location of the shear band depends 
on the details of the initial fluctuations in x and therefore 
seems arbitrary on macroscopic scales. This is similar to 
what is seen in Shi, et al. for simulations with periodic 
boundary conditions. However, in other MD simulations 
and in experiments on bulk metallic glasses the bound- 
ary conditions are non-trivial and most likely play an 
important role in nucleating shear bands. For example, 
in simulations by Varnik et al. [l9[ the system is driven 
at the boundary by a rough, rigid layer and the material 
consistently develops shear bands at the boundary. Ad- 
ditional simulations by Falk and Shi [StJ show that shear 
bands are more likely to occur in the presence of surface 
defects. 

Boundary conditions may also influence the long-time 
behavior of shear bands. The long-lived, localized so- 
lutions to the STZ equations are not truly stationary 
states - at long times the shear band will diffuse to 
cover the entire strip and the system will become ho- 
mogeneous. Figured] is a plot of the steady state stress 
vs. the natural log of the imposed strain-rate for STZ 
solutions. This function is not multivalued, indicating 
that the only steady state solution to the STZ equations 
with periodic boundary conditions is the spatially ho- 
mogeneous solution x(y) = Xoo, s — Sf. Even exclud- 
ing diffusion, the STZ equations will evolve to a homo- 
geneous solution because the plastic dissipation rate is 
positive everywhere, which means the effective tempera- 
ture is everywhere driven towards Xoo, albeit at different 
rates. However, if the effective temperature field is speci- 
fied at the boundaries (i.e, the boundary causes ordering 
and fixes the disorder temperature there) then it is pos- 
sible for solutions to have inhomogeneous steady states. 
This is similar to results for simulations with rough, rigid 
boundaries where the material develops stationary local- 
ized states [19j. More work is needed to understand the 
relationship between shear band development and mate- 
rial boundaries, and STZ theory should provide an excel- 
lent framework for these investigations. 

The athermal STZ description provides a simple start- 
ing point for understanding the dynamics of shear band 
formation. As discussed in the previous section, a ther- 
mal description is likely a better model for the MD sim- 
ulations at different strain rates. Adding thermal relax- 



17 




Log 10 (Strain rate/T ) 



FIG. 8: STZ solution for the steady state stress vs. Logio 
of the imposed total strain rate. The curve is single- valued 
throughout, indicating that the system does not posses steady 
state coexisting regions with two different phases unless im- 
posed by the boundary conditions. 

ation to the STZ equations generates a picture of local- 
ization similar to those found in complex fluids, where 
the stress vs. strain rate curve is multi- valued [481 ]. As 
discussed by Langer [2|| and Foglia [4l[ , relaxation terms 
in the equation for the effective temperature generate so- 
lutions where the steady state stress is a multi-valued 
function of the strain-rate. In steady state, two regions 
with different strain rates can coexist: one region which 
is nearly jammed except for thermal rearrangements and 



a second region where STZs are constantly created and 
flipped, resulting in plastic flow. 

This is different from the athermal theory where all 
of the material deforms due to shear- induced transitions, 
though the transitions occur at different rates. While 
athermal analysis provides a simple framework for under- 
standing mechanisms leading to shear band nucleation, 
an STZ formulation with thermal relaxation may further 
enhance our understanding of steady-state dynamics and 
connect mechanisms for localization in amorphous solids 
to those in complex fluids. 
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